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Abstract 

Antiferromagnetic Ising spins on the scale- free Barabasi- Albert network are studied via the Monte 
Carlo method. Using the replica exchange algorithm, we calculate the temperature dependence 
of various physical quantities of interest including the overlap and the Binder parameters. We 
observe a transition between a paramagnetic phase and a spin glass phase and estimate the critical 
temperature for the phase transition to be T ~ 4.0(1) in units of J/kB, where J is the coupling 
strength between spins and ks is the Boltzmann constant. Using the scaling behaviour of the 
Binder parameter, we estimate the scaling exponent to be v ~ 1.10(2). 
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I. INTRODUCTION 




In the last few years, the study of complex networks has found relevance in various fields 
including sociology, ecology, biology, economics and physics. In these networks, vertices do 
not have homogeneous links or connectivities. A particularly relevant structure found in 
several empirical studies is the so-called scale-free network (SFN), which is characterized by 
the power law distribution of the degree of connectivities, P(k) ~ fc~ 7 , with k the number 
of links for a node, and 7 the decay exponent of the distribution. A network with 7-^0 
has nodes with a relatively homogeneous number of links (somewhat resembling the case 
on regular lattices), while large 7 corresponds to the existence of "very famous" nodes (or 
hubs), i.e., those having direct links to the majority of vertices. 

Many networks realized in Nature show scale-free structure. Some examples studied 
include food webs jlj], power grids and neural networks Q, [|, cellular networks 
contacts Hi, Internet routers |Q|, the World Wide Web l7f, actor collaborations 

n n 

the citation network of scientists 10J and the stock market 

In addition to the scale-free behaviour, these networks are characterized by a high clus- 
tering coefficient, C, in comparison with random graphs [12]. The clustering coefficient, C, 
is computed as the average of local clustering, Cj, for the i th node, defined as 

Ci = (1) 

where Z{ is the total number of nodes linked to the site % and yi is the total number of 
links between those nodes. As a consequence both and C lie in the interval [0,1]. The 
high level of clustering found supports the idea that a herding phenomenon is a common 
feature in social and biological communities. The parameter C also represents the density 
of triangles, that is of elementary cells, associated with the network. 

Numerical studies on SFNs have demonstrated how topology plays a fundamental role in 
infection spreading fl^ . opinion formation in large communities [b^ and tolerance against 
random and preferential node removal 14 , 3| • A detailed description of the pro gres s in thi s 



emerging field of statistical mechanics can be found in the recent reviews ofrMs. j i4li2l,Q • 
The aforementioned empirical findings have inspired physicists to investigate the dynam- 
ics of standard models in the new case where the interactions between elements are described 
by complex networks. These include the study of various magnetic models such as the Ising 
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model. An intriguing issue concerns how the unusual topology acts to influence the cooper- 



ative behaviour of the spins. Studies of t 



several theoretical techniques 19, 
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re ferromagnetic (FM) Ising model on a SFN, usin 
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22], 



, 2^ including the Monte Carlo (MC) method 
have found the robustness of ferromagnetic ordering against thermal fluctuations for the 
degree distribution exponent 7 < 3. This result is actually intuitive if we notice that, as 7 



gets smaller, nodes at the edge of the network will generally have more connections. In this 
situation, the system resembles the FM Ising model on a regular lattice which exceeds the 
lower critical spatial dimension, di = 2. There the ordered phase is very robust against ther- 
mal fluctuations. However, for the antiferromagnetic (AF) case with a SFN, the situation is 
different. 

Two factors come to play a central role in the dynamics of the AF-SFN model; namely 
the competition induced by the AF interaction in the elementary triangles of the network 
and the randomness related to the non-regular connections. The abundance of elementary 
triangles in the network leads to frustration, as, for example, only two of the three spins can 
be anti-aligned. More generally, frustration refers to the inability of the system to remain 
in a single lowest energy state (ground state). These ingredients lead the AF SFN to belong 
to a class of randomly frustrated systems commonly referred to as spin glasses (SGs). 

Most studies of SGs have been performed on regular lattices. These studies have shown 
that frustration and randomness are the key ingredients for SG behavior, characterized by 
a frozen random spin orientation at low temperatures j^ . Spin glasses on a SFN with 
mixed AF and FM bonds have been investigated recently by Kim et al. j^j]. They found, 
for 7 < 3 and even distributions of the two kinds of interaction, that the system is always 
in a SG state for any finite temperature. A study of the pure AF Ising model on a SFN is 
of great theoretical interest since, despite the homogeneity of the bonds, it inherits all the 
characteristics of a SG from the random frustration related to its geometry. General reviews 
on SG systems can be found in Refs. j^J. 

In this paper we consider the AF Ising model on a SFN, more precisely the Barabasi- 
Albert (BA) network with tunable clustering . Using the replica exchange algorithm 
of the Monte Carlo method, we calculate the order parameters of spin glass behaviour, the so- 
called overlap parameter and its distribution. For an accurate determination of the critical 
temperature, we also evaluate the Binder parameter. The paper is organized as follows: 
Section HH describes the model and the method. The results are discussed in Section IIHI 
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FIG. 1: (Color online). Example of a scale- free network. The number of nodes is 500 with clustering 
probability 9 = 0.9 and mo = m = 2. The number of nodes has been kept small in order to preserve 
the clarity of the plot. Note that, for such small networks, a large scale invariant range is obtained 
only if one considers the ensemble average over several realizations. This plot has been realized 
with the Pajek software [2^ . 

Section IIVI is devoted to the concluding remarks. 

II. MODEL AND SIMULATION METHOD 
A. The model 

In order to create the scale-free network topology we make use of the Barabasi-Albert 
model P|. This is based on two main considerations: (i) linear growth and (ii) preferential 
attachment. In practice the network is initialized with m disconnected nodes. At each step 
a new node with m edges is added to the pre-existing network. The probability that an edge 
of the new node is linked with the ith node is expressed by Tl(ki) = kij Y^j kj. The iteration 
of this preferential growing process yields a scale free network, where the probability of 
having a node with k connections is P(k) ~ fc -7 with 7 = 3. This is an interesting value. 
In the thermodynamic limit, the second moment of the distribution diverges, (k 2 ) = 00, for 
7 < 3. This leads to peculiar properties of theoretical models in this range of 7 values (l8|. 
In the present work we focus on the case in which 7 = 3 and the divergence of (k 2 ) is 
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logarithmic. An extensive investigation of the phase space for the AF model on SFN is left 
for future work. 

It is also worth noting that the Barabasi- Albert model cannot reproduce a high clustering 
coefficient. In fact, the value of this coefficient depends on the total number of nodes, N, in 
the network and in the thermodynamic limit, N — > oo, C — > 0. 

In the AF Ising system the average cluster coefficient, C, plays a fundamental role in the 
dynamics. In fact, it represents the average number of triangles per node and, as a result, it 
is directly related to the degree of frustration in the network. In order to keep this parameter 
constant, on average, with the size of the network, we introduce a further step in the growth 
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In this case, if the 



process, namely the triad formation proposed by Holme and Kim 
new added node is linked with an older node, i, having other links, then with a certain 
probability, 9, the next link of the new node, if any remain, will be added to a randomly 
selected neighbour of node i. This method of introducing friends to friends, while preserving 
the scale-free nature of the networks with 7 ~ 3, generates high clustering coefficients that 
do not depend on N. The only tunable parameter that changes the value of the clustering 
coefficient is the clustering probability 9. An example of a SF network generated with this 
algorithm is shown in Fig. ^ for 500 nodes. 

We simulate various sizes of the network with many different realizations and investigate 
the scaling behaviour of the various physical quantities we are interested in. All the simula- 
tions have been carried out fixing 9 = 0.9, corresponding to an average clustering coefficient 
of C ~ 0.39, close to the value found in many real systems On each SFN constructed 
at the beginning of the simulation, we assign to each vertex an Ising spin, and to each link 
an AF interaction. The Hamiltonian can be written as follows 

H = - J n s i s i ■ ( 2 ) 

(«> 

Here the summation is performed over the connected spins Sj and Sj occupying sites % and 
j, respectively. The coupling interaction J^- = J = — 1 is AF. As previously mentioned, 
each vertex with the local cluster coefficient Cj > together with its neighbours, compose 
elementary triangles. Due to the AF interactions the local system is frustrated. 

It is worth pointing out that C is related to the degree of frustration of each network. Due 
to the probabilistic algorithm used for their construction, the value of C fluctuates around 
a mean value from one network to the next and, therefore, provides a source of randomness 



that, as we will see, gives rise to the spin glass properties of the model. This probabilistic 
growth is not shared by other algorithms which use recursion formulas to generate scale-free 
structures, such as, for example, the Apollonian networks 28[. In this case, once one fixes 
the number of iterations of the algorithm, which is proportional to the number of nodes of 
the final network, one also fixes its topology. The element of randomness is therefore missing 
in the Apollonian procedure. 

As a random system, each realization of a network of size iV will differ in the "structure" 
of connectivities. Therefore, in order to have reliable statistics, we average over many re- 
alizations of the SF network for each specified size. The system sizes that we simulate are 
N = 1024, 2048, 4096, and 8192. In general, one takes into account more realizations for 
small system sizes and less for large system sizes as the latter tend to self-average. However, 
since the self-averaging of physical quantities for larger system sizes is interfered with by 
the increase of ground state degeneracy, we do not take less realizations. Instead all phys- 
ical quantities of interest for each system size are averaged over 1000 network realizations. 
Moreover, for each realization of the network, we fix mo = m = 5, corresponding to a coor- 
dination number on a regular lattice of approximately 10. In the thermodynamic limit, the 
average connectivity for the BA network is (k) = 2m = 10, emphasizing the fact that we 
are implicitly dealing with a high dimensional system. 

Another peculiarity of SF networks is the existence of a broad distribution of "hubs", 
that is nodes with a large number of connections, k. The energy difference in a spin flip 
actually depends on the number of connections of the spin itself, AEi = -2siYl}Usy Thus 
in the AF case for the ith spin with k{ connections, the hubs are more likely to "freeze" 
into a particular configuration compared to the nodes with just a few links. This property 
resembles the spin glass behaviour of particular alloys where some elements freeze into a 
particular orientation at a higher temperature than others. 



B. Simulation method 



The calculation of the thermal averages of the physical quantities of interest is performed 
using the replica exchange MC method [26] . In this method the evolution of M replicas, each 
in equilibrium with a heat bath of inverse temperature j3 m for the m th replica, is simulated 
in parallel. Given a set of inverse temperatures, {/?}, the probability distribution of finding 
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the whole system in a state {X} = {X±, X 2 , . . . , Xm} is 

M 

P({X,f3})= J] P(X m ,Pm), (3) 

m=l 

with 

P(X m ,/3 m ) = Z(/? m )- 1 exp(-/? m i/(X m )), (4) 

and Z(/3 m ) is the partition function at the m th temperature. We can then define an exchange 
matrix between the replicas in our Markov chain, W(X m , /3 m \X n , /3 n ), that is the probability 
to switch the configuration X m at the temperature (3 m with the configuration X n at (3 n . By 
using the detailed balance condition, required to keep the entire system at equilibrium, on 
the transition matrix 

P{- ■ ■ , \X m , P m }, . . . , {X n , p n }, . . .) • W(X m , (3 m \X n , p n ) 
= P{. . . , {X n , /? m }, . . . , {X m , p n }, ...)■ W(X n} f3 m \X m , p n ), (5) 

along with Eq. (J3J), we have that 

Y R I Y R \ 

exp(-A), (6) 

where A = (/3„ — f3 m )(H(X m ) — H(X n )). With the above constrains we can choose the 
matrix coefficients according to the standard Metropolis method and, therefore, we have 

x f 1 if A < 0, 

W(X m ,(3 m \X n ,{3 n ) = { (7) 

[ exp(-A) if A > 0. 

In our simulation we restrict the exchange to temperatures next to each other; that is, we 
consider only the terms W(X m , /3 m \X m+ i, f3 m +i). This choice is motivated by the fact that 
the acceptance ratio decays exponentially with (/3„ — p m ). 

The replica exchange method is extremely efficient for simulating systems such as spin 
glasses, that can otherwise become frozen in some particular configuration at low temper- 
atures when using a standard Metropolis algorithm for the configuration update. In this 
case, as we lower the temperature, the system can become trapped into a local minimum of 
the free-energy where the barriers are so high that the time required for the system to move 
to another allowed region of the configuration space diverges to infinity as a function of the 
system size. If the system is trapped in a local minimum then the ergodicity condition is not 
fulfilled anymore and the measure that one makes become biased by the particular region 



W(X m ,(3 m 


\X n , P n ) 


w(x n ,p m \ 


Xfm Pn) 



of the configuration space that is being sampled. By using the exchange replica method, 
instead, we keep switching the temperatures between the M copies of the system and, as 
long as the higher temperature is in a hot phase (where, the system can easily explore 
all the configuration space), then we are in principle able to explore all the configuration 
space also for the lower temperatures. Another advantage of this method is that the replica 
exchange reduces drastically the temporal correlation in the system dynamics at each tem- 
perature. This enables one to collect more independent measures for the thermal averages 
of the physical quantities and, therefore, reduces the uncertainty. 

It is important to stress that, before starting the actual simulations, some care is required 
in selecting the set of inverse temperatures, {/?}. In fact, the method is efficient only when 
a fairly large transition probability is maintained in the range of interest. From Eq. we 
can see that, in the hot phase, temperatures can be more coarsely spaced while in the cold 
phase the temperatures need to be closer to each other. An optimal set of temperatures can 
be obtained by iterating, in preliminary runs, the following map j^fil ]: 

A = A " (8) 

Pm = P m -1 + (Pm ~ Pm-l) ' Pm/c, 

where p m is the acceptance ratio for the switch between two configurations at the mth 
temperature and c = J2m=iPm/ (M — 1) is a normalization factor. The initial value for the 
set {/5} is uniform in the interval of interest and we ensure that /3% belongs to the hot phase. 
For each iteration of the map, a run of a few thousand MC steps is carried out to calculate 
the acceptance ratios, p m , which are then plugged into Eq. (JHJ) in order to obtain a new 
set of inverse temperatures. After a few iterations, the map of Eq. (jHJ) converges to a fixed 
point, {/?*}, which sets the values of the temperatures to be used in our simulations. 

In using this method, we define a "local" MC (LMC) update as a MC update for each 
spin of each replica, either consecutively through all elements of the network or randomly. 
Given that we can group the inverse temperatures in even and odd pairs, (J3 m ,ft m+ x), after 
each LMC update we alternate attempts to switch configurations from one temperature to 
the next. According to this procedure, we define a Monte Carlo step (MCS) as a LMC plus 
a half (m odd or even) exchange trial. 

For each realization of the network we start from a random configuration of the spins and 
then perform 10 3 LMC updates in order to reach thermal equilibrium. After this transient 
period, we run the simulation for 3 x 10 5 MCSs while taking a total of 6 x 10 4 measures for 
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the thermal averages, that is one every 5 MCSs (temporal correlations are lost very quickly 
by using the replica exchange method). We consider low temperatures in a search for the 
possible existence of a phase transition. The thermal averages obtained for each network 
are then averaged over the ensemble of networks. In the following, we indicate (...) as the 
thermal average and [...].„ as the ensemble average. The statistical errors in the plots, where 
reported, are calculated via the bootstrap method. 



III. RESULTS AND DISCUSSION 

A. Spatial correlations and specific heat 

As a first step we investigate the extent of spatial correlation of the spins in the SF 
network by making use of the spatial autocorrelation function which is defined on a regular 
lattice as 



i 

av 



(9) 



where Ld is the total number of pairs at distance r and depends just on the dimension 
considered. In a SF network the situation is more complicated since there may be several 
paths leading from a certain node to another. We then define r as the minimum path 
between two nodes and the denominator of the Eq. (JOJ) becomes dependent on r. The 
results, averaged over 50 configurations, between the temperatures of T = 5.0 and T = 2.1 
are shown in Fig. El for iV = 1024. All the temperatures in the present paper are expressed 
in units of J/ kg, where J is the coupling strength between spins and ks is the Boltzmann 
constant. 



In order to give a better interpretation of the plot in Fig.Elwe remind the reader about an 
important propriety of SF networks; that is their "small world structure" . The "hubs" , in 
fact, play a fundamental role in linking sites otherwise very distant. Moreover, the average 



path length increases just logarithmically with the size of the network [16l 1 171 ]. In the plot 
of Fig. |2l for N = 1024 nodes, an upper limit of r = 6 is encountered. While all the 50 
configurations reach r = 6, only a few networks exceed this limit. 

The plot emphasizes how neighboring spins, on average, tend to be anti-correlated, as 
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FIG. 2: (Color online). Spatial autocorrelation, £(r), for N = 1024 averaged over 50 network 
configurations for temperatures between T = 5.0 and T = 2.1. The plot shows that next neighbour 
spins tend to be anti-parallel as in standard AF Ising model. The AF interaction in the triangular 
units of the system results in high frustration. Note that the number of nodes at large distances is 
much smaller than the ones at smaller distances and so the average calculated for r = 5 and r = 6 
includes just few samples. This is a consequence of the "small-world" effects in SF networks. 




FIG. 3: Specific heat, C„, as a function of the temperature and system size. The plot has been 
obtained by averaging over 50 network configurations for each N . Note that the specific heat does 
not scale with the size of the system. 
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expected in the AF case. The autocorrelation decreases with the distance from the node 
under consideration. The temperature dependence is also in accord with the expectations. 
The absolute value of the correlation decreases with increasing temperature and vice versa. 
Indeed, the highest and lowest temperatures form a perfect boundary for all the curves. This 
is an expected result, since thermal effects always tend to reduce the correlation between 
the spin interactions. 

We also study the behaviour of the specific heat, C u , defined as follows 



where is the Boltzmann constant. Although no singularity is expected for this quantity 
in the spin-glass transition, it is interesting to compare its behaviour with other studies. 
The dependence of the specific heat on temperature is reported in Fig. El The statistical 
errors, in this smaller than the size of the symbols and therefore are not reported. A 

common Schottky peak of the specific heat for a finite system is observed at the temperature 
of T ~ 2.0 independent of the system size. Below this point, we found that C v decreases 
and goes to zero as T — > 0. 

This behaviour follows from simple entropy considerations. In fact, since we are dealing 
with a finite Ising system, the entropy is bounded at each finite temperature as well, 



and, necessarily, C v — > for T — > 0. 

The next section is dedicated to study of the SG behaviour and the phase transition of 
the system. In order to achieve this task, we evaluate the corresponding order parameters, 
the overlap parameter and the Binder parameter. 

B. Observing spin glass behaviour 

With the presence of frustration and randomness in the AF-SFN model, we expect to 
observe a spin glass transition, i.e., a transition from a temporal disordered to a temporal 
ordered phase at low temperatures. 

This feature is not shared by the so-called fully frustrated systems This type of 

transition might be characterized by the order parameter such as that suggested by Edward 




(10) 




(11) 
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and Anderson 301, defined as follows 

" f 



qEA 



(12) 



However, an ergodic Markov chain of a system having Z 2 symmetry will ensure the thermal 
average of the ith spin vanishes. Therefore a finite value of this measure simply reflects the 
non-ergodicity in the MC update. 

A more appropriate quantity that is often used to characterize the SG state is the overlap 
parameter, q, defined as j^H Q 

1 = ^4 a) sf\ (13) 

i 

where the superscripts a and (3 denote two copies of the same configuration of connectivity 
at the same temperature. The actual value of q is extracted from both the thermal and 
disorder average, [(...}] av . 

Using the replica exchange MC simulation, the two copies, a and /3, are allocated at 
each temperature of the parallel tempering. This means, if the measurement is performed 
on M points of temperatures, there are M pairs of replicas. The Metropolis spin update 
is performed on each node for every MC step. As a part of the equilibration steps of the 
algorithm described in the previous section, we exchange two a (and (3) replicas of neigh- 
boring temperatures, according to a certain probability. Then, for each temperature, the a 
and (3 replicas are superimposed every 5 MCSs in order to measure the overlap parameters, 
as defined in Eq. fT3|) . 

In particular, for the Ising system, due to the Z 2 symmetry, it is important to evaluate 
the absolute value of the order parameter, 

1 



(14) 



to overcome the implication of the Z 2 symmetry of the Hamiltonian, that is the config- 
urations Si and — Si have equal Boltzmann weights. That is, if the system is at thermal 
equilibrium and if we take quite long MCS then the usual q should average to zero. The 
existence of a spin glass phase is indicated by the convergence of \q\ to a finite value as we in- 
crease the network size. At the same time, a convergence of |g| to zero at high temperatures 
is anticipated. In the latter case the system is in the paramagnetic phase. 
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FIG. 4: Temperature dependence of the overlap parameter, q, for different system sizes N. The 
increasing value of q at low temperatures indicates a SG phase. For a given network size, 1000 
realizations of the SFN are averaged over. 

The temperature dependence of \q\, resulting from the simulations, is shown in Fig.EJ The 
existence of a SG phase is indicated by the finite value of |g| in the low temperature region, 
and the approach of \q\ to zero at higher temperatures associated with the paramagnetic 
phase. For high temperatures and large networks, \q\ is approaching zero in accord with the 
thermodynamic limit where \q\ = |33j. 

The existence of these two different phases can also be observed from the distribution of 
q, as shown in Fig. El For higher temperatures we observe simple Brownian fluctuations of 
the values of q, leading to a singly peaked Gaussian distribution characteristic of a param- 
agnetic state. By decreasing the temperature, the distribution spreads out, reflecting the 
increasing number of metastable disordered states associated with a substantial frustration. 
At lower temperatures the distribution develops double peaks reflecting the Z 2 symmetry 
and a finite value of \q\, representative of the SG phase. We note that the shape of the 
observed distribution at low temperatures is different from that of the conventional Ising 
system where the double peaks approach delta-like double peaks reflecting a simple doubly 
degenerate ground state |34j ]. 

An accurate evaluation of critical temperature of the phase transition is achieved via the 




13 



Binder parameter defined as follows 

Here (q 2 ) and (q 4 ) are respectively the second and the fourth cumulant moment of q. In 
this calculation, in order to avoid systematic correlation errors that could bias the results if 
we were evaluating this average over g L directly 3^|, the second and fourth order cumulants 
are averaged prior to taking their ratio. The Binder parameter is constrained in the range 
< < 1- At high temperature, where thermal fluctuations overcome all cooperative 
interaction, the system is expected to exist in the paramagnetic phase where there is no 
spatial autocorrelation. As a result, the distribution of q should be Gaussian centered at 
q = 0. In this case the ratio of the cumulants, (g 4 )/ (q 2 ) 2 — > 3, resulting in g L — ► 0. At low 
temperatures, the cooperative interaction becomes dominant and the ratio of the cumulants 
approaches unity so that gL — > 1. 

Fig. |U1 (inset) displays the temperature dependence of the Binder parameter for a variety 
of network sizes. A spin-glass state is observed for lower temperatures where the Binder 
parameter deviates from zero, and increases with the system size while approaching to 1. 
In the thermodynamic limit, we expect gi — > 1 just below the critical temperature. A 
crossing point in the size dependence of g^ indicates that the critical temperature for the 
SG phase transition is T ~ 4.0. Fig. |U] indicates that for temperatures above T ~ 4.0 the 
Binder parameter, while remaining always above zero, does indeed order in an opposite 
manner indicative of a genuine crossing of the curves and in accord with a genuine spin 
glass transition at finite temperature. This feature which is not observed for uniformly 

n 

distributed AF and FM bonds, as T c = oo in the thermodynamic limit [24|. However, the 
value of the transition temperature is not determined with high accuracy by the crossing of 
the Binder parameter. In fact, finite size effects seem to slightly distort the tendency for 
very small networks, as in the case of N = 1024. At the same time, the statistical errors 
in the paramagnetic phase for large networks, see N = 8192, appear to be significant and 
some points are scattered. 



A more accurate estimate of the critical temperature, T c , for finite size systems can be 
obtained using scaling arguments. For a SG system, the Binder parameter depends on the 
system size L as 
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FIG. 5: (Color online). The distribution of q at various temperatures for different system sizes, 
including (a) N = 1024, (b) N = 2048, (c) N = 4096 and (d) N = 8192. 
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FIG. 6: Scaling behaviour of the Binder cumulant, g^, for different system sizes. Each system size 
is averaged over 1000 realizations of the network configuration. 
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FIG. 7: Scaling plot of the data illustrated in Fig. El fitted to Eq. 1161 



g L = ~g L [{T-T c )L 1 /% (16) 



cumulant does not depend on L. 
the number of nodes N Q, 12, 3, 3 



being z/ > the spin glass correlation length exponent, implying that at T c the Binder 

7 or the SFN, the system size scales logarithmically with 
and therefore we take L = log(iV). This slow increase 
in the diameter of the system, as well as the average path length, is a manifestation of the 
"small- world" property of this network, induced by the presence of a large number of highly 
connected hubs which create shortcuts between the nodes. An important implication of 
this feature is that we cannot embed the network in any finite dimensional lattice: we are 
implicitly dealing with a high dimensional system. The correlation length, in this case, is 
still well defined although its value gets close to the densely-connected, mean field limit as 
we increase the average connectivity of the nodes, (k) = 2m. 

The parameters T c and v are determined by constraining the temperature dependence of 
the Binder parameter for each network size to lie on a single curve. The curves following 
the scaling bahaviour of Eq. ()16|) are shown in Fig. [7| From this fit we estimate the critical 
temperature T c ~ 4.0(1) and the exponent of the SG correlation length v ~ 1.10(2). It is 
important to underline that this kind of behaviour is not observed for an AF system on a 
regular triangular lattice. 
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IV. CONCLUDING REMARKS 



In summary, we have investigated the antiferromagnetic Ising model on a Barabasi-Albert 
scale-free network using the replica exchange Monte Carlo method. Through the calculation 
of the overlap parameter we observe spin glass behaviour at low temperatures. Using the 
scaling behaviour of the Binder parameter the critical temperature separating the SG and 
the the paramagnetic phases is found to be T c = 4.0(2) with a scaling exponent of SG 
correlation length v ~ 1.10(2). Such behaviour is not observed for the AF Ising model on 
regular triangular lattices. Hence the topology of the interactions plays a critical role in the 
dynamics of the system. 
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